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ABSTRACT 

The dynamics of the partially ionised solar photosphere-chromosphere can be de- 
scribed by a set of equations that are structurally similar to the magnetohydrody- 
namic equations, except now the magnetic field is no longer frozen in the fluid but 
slips through it due to non-ideal magnetohydrodynamic effects which manifests itself 
as Ohm, ambipolar and Hall diffusion. Macroscopic motion of the gas is widespread 
throughout the solar atmosphere and shearing motions couple to the non-ideal effects, 
altogether destabilising low frequency fluctuations in the medium. The origin of such 
non-ideal magnetohydrodynamic instability lies in the collisional coupling of the neu- 
tral particles to the magnetized plasma in the presence of sheared background flows. 
Expectedly the maximum growth rate and most unstable wavenumber depend on the 
flow gradient and ambient diffusivities. 

The orientation of the magnetic field, velocity shear and perturbation wavevector 
play crucial role in assisting the instability. In the presence of only vertical field and 
vertical wavevector, ambipolar and Ohm diffusion can be combined together as Ped- 
ersen diffusion and causes only damping; in this case only Hall drift in tandem with 
shear flow drives the instability. However, for non- vertical field and oblique wavevector, 
both ambipolar diffusion and Hall drift assist the instability. 

We investigate the stability of a magnetic element in the network and internetwork. 
The shear scale is not yet observationally determined and thus, assuming typical shear 
flow gradient ~ 0.1 s _1 we show that the magnetic diffusion shear instability grows 
over a minute. Thus, it is plausible that network-internetwork magnetic elements are 
subject to this fast growing, diffusive shear instability. This fast-growing instability 
could play important role in driving low frequency turbulence in the photosphere- 
chromosphere plasma. 
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1 INTRODUCTION 

The solar atmosphere is host to a number of highly ener- 
getic events such as coronal mass ejections (CMEs), flares, 
prominences, coronal heating. The huge reservoir of energy 
in the solar atmosphere is related to the magnetic field orig- 
inating in the convection zone where shearing motion of 
field line footpoints can stretch and twist the anchored field 
(Parker 1979). A variety of magnetohydrodynamic (MHD) 
waves are generated in the corona due to convective mo- 
tion in the photosphere. The convective shearing motion not 
only produces Alfven , fast or, slow magnetoacoustic waves, 
but, also, brings topologically disparate parts of magnetic 
configuration closer together resulting in the formation of 
current sheets. All in all, a tiny fraction of convective en- 
ergy carried by the waves to higher altitude may suffice to 



heat the coronal plasma to high temperatures. This simple 
physical picture of wave excitation, propagation and ensu- 
ing coronal heating is quite attractive since it ties the heat 
transport in the solar corona and heliosphere to the ultimate 
source of energy - the shearing photospheric convective mo- 
tions. Therefore the investigation of wave propagation and 
concomitant heating of coronal plasma in the framework of 
MHD has been very popular topic in solar physics (Priest 
1987; Parhi et al. 1997a,b, 1998; Hasan & Kalkofen 1999; 
Goedbloed & Poedts 2004; Aschwanden 2009). 

Most solar atmospheric heating, with possible excep- 
tion of flares, takes place in the chromosphere. The chromo- 
sphere is highly stratified and extends up to about nine pres- 
sure scale heights, i.e. about 2 Mm above the photospheric 
surface. The lower chromosphere is threaded by strong (~ 
kiloGauss) vertical flux tubes located in the network regions 
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where they are observed as bright points. These tubes, which 
have a low filling factor (< 1%) near the footpoint in the 
photosphere, expand to fill about 15% of the lower chromo- 
sphere (~ IMm, where Call emission features are observed 
in H and K lines) before filling the entire atmosphere forming 
a canopy in the chromosphere. The quite solar internetwork 
region is also magnetised, with patches of kG field concen- 
tration (de Wijn et al. 2009; Almeida et al. 2010) and an 
order of magnitude smaller field everywhere else. Observa- 
tion suggests that localised active regions that emit ~ 80% 
of the coronal radiative loss near solar maximum contain 
plasma of chromospheric origin (Aschwanden 2001). This 
fact raises the possibility that the same mechanism that 
transports mechanical energy from the convection zone to 
the chromosphere to sustain its heating rate also supplies 
the energy needed to heat the corona and accelerate the 
solar wind. 

The absence of a unified theoretical framework to ad- 
dress the heating problem of the chromosphere-corona is 
partially related to the overuse of ideal MHD, from the deep 
interior of the sun to the corona. This approach is valid 
only for the fully ionised plasma in which field and fluid 
are well coupled to each other, a far cry from the low tem- 
perature photosphere-chromosphere where non-ideal MHD 
effects such as Hall, ambipolar and Ohm diffusion domi- 
nate (Pandey & Wardle 2006, 2008; Pandey et al. 2008; 
Khomenko & Collados 2012; Sykora et al. 2012; Pandey & 
Wardle 2012a). It is well known [Vernazza et al. (1981); here- 
after VAL81] that the number of plasma particles in the par- 
tially ionised solar atmosphere, particularly below 2.5 Mm 
is small compared to the neutrals and thus, the plasma fluid 
is not frozen in the magnetic field but can slip through it 
due to frequent collisions with the neutrals (Mestel & Spizer 
1956). 

The inclusion of neutral dynamics not only destroys 
the economy and simplicity of single fluid MHD descrip- 
tion of fully ionised plasma but very concept of well de- 
fined flux tube may be difficult to define in the multi-fluid 
framework. Furthermore, high frequency Alfvenwave may 
not survive the collision dominated photosphere and chro- 
mosphere (Goodman 2004; Leake et al. 2005; Vranjes et al. 
2007; Arber et al. 2007; Vranjes et al. 2008; Goodman 2011). 
This could have been anticipated on the ground that in the 
photosphere (< 500 km) and chromosphere (< 2500 km) the 
plasma number density is much smaller than the neutral 
number density and thus, high frequency (with respect to 
the ion-neutral collision frequency) MHD waves are severely 
damped by collisional dissipation of the wave energy. A way 
out of this difficulty is to retain the MHD momentum equa- 
tion and include collisional effects by modifying the induc- 
tion and energy equations via conductivity tensor (Erdelyi 
& James 2004; Leake & Arber 2006), an approach often em- 
ployed to study the lower ionosphere of the Earth . However, 
derivation of the time-independent conductivity tensor re- 
quires the neglect of time derivatives in the electron and ion 
momentum equations, i.e. d e ^i/dt ~ <C uj C e,ci, implying 
plasma dynamical response frequency, uj e< i is much smaller 
than their respective gyro-frequencies, oj C e,ci. Therefore, ne- 
glecting plasma inertia, a linear relationship between the 
electric field E and plasma current J can be easily derived 
E = er ■ J where <x is the time-independent conductivity 
tensor. However, equation of motion in the MHD frame- 



work assumes w< ~ oj ci . In fact ion carries the inertia of 
the fluid. Therefore, on the one hand, the time-independent 
conductivity tensor implies u) et i <C <x) C e,ci, on the other hand, 
MHD momentum equation requires u)i ~ uj c i. Clearly, inves- 
tigation of the collisional effects by merely modifying the 
induction and energy equation in the MHD framework is 
highly unsatisfactory. 

Any realistic model of the solar atmosphere must re- 
flect two basic observational facts: (a) the magnetic field 
distribution on the solar surface is not continuous but is or- 
ganised in the network and internetwork elements. Whereas 
network field (> kG) is predominantly vertical organised 
in the flux tubes (diameter < 100 km) located in the in- 
tergranular lanes, the internetwork field [G (few G — kG)] in 
the interior of supergranule cells is mostly horizontal (Hasan 
2009; Lites et al. 2008) 1 ; (b) plasma in the photosphere- 
chromosphere is weakly ionised [with fractional ionisation, 
which is a ratio of the electron and neutral number densities 
X e = n e /n n ~ 10~ 4 , VAL81]. 

Both active and quite phases of solar atmosphere are 
highly dynamic and consist of convectively driven vortices 
and flows of various spatial and temporal scales (Bonet 
et al. 2008; Wedemeyer-Bohm & Voort 2009; Balmaceda 
et al. 2010; Bonet et al. 2010). Most vortices are small 
(< 0.5 Mm) with average size ~ 241 ± 25 km and typi- 
cal lifetime ~ 3 — 5 min. though large vortices ~ 20 Mm 
with lifetime > 20 min. have also been observed (Attie et al. 
2009) . The bright points associated with the vortex motion 
in the intergranular lane moves with typical speed < 2km/s 
(Wedemeyer-Bohm & Voort 2009). The magnetic fields ap- 
pear to play crucial role in mediating vortex motion in the 
photosphere and chromosphere (Steiner & Rezaei 2012). 

Rotation has been invoked in the past to explain the sta- 
bility of flux tubes (Schussler 1984). Models of spicules also 
use concept of rotating flux tubes (Kudoh & Shibata 1997). 
The numerical simulation of solar convection displays turbu- 
lent vortex flows at intergranular lanes (Zirker 1993; Stein 
& Nordlund 1998). The vorticity generation near bound- 
aries of granules have also been seen in the numerical sim- 
ulation of photosphere (Nordlund et al. 2009; Muthsam et 
al. 2010). The formation of small-scale, intergranular vor- 
tices suggests that vorticity is formed due to interaction of 
the photospheric plasma with ambient magnetic field in the 
intergranular lanes (Moll et al. 2011; Shelyag et al. 2011). 
The high resolution simulations of non-ideal MHD effects 
shows that Hall effect generates out-of-plane velocity fields 
with maximum speed ~ 0.1 km/s at interface layers between 
weakly magnetized light bridges and neighbouring strong 
field umbral regions (Cheung & Cameron 2012). To sum- 
marise, both observation and numerical simulation points 
to the presence of shear flow at various spatial scales in the 
solar photosphere. 

The large scale shear flow acts as a source of free energy 
in the solar plasma that can easily destabilise waves. For 
example, shear driven Kelvin- Helmholtz instability (KHI), 
which converts shear flow energy into vortex kinetic en- 
ergy, is invoked to explain the instability of the flux tubes 
(Kolesnikov et al. 2004; Soler at al. 2010). Although solar 



1 Internetwork field could be as well isotropic (Almeida and 
Gonzalez 2011), but see also Steiner & Rezaei (2012). 
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atmosphere may be susceptible to KHI (Karpen et al. 1993; 
Kolesnikov et al. 2004; Soler at al. 2010; Zaqarashvili at 
al. 2010; Ofman & Thompson 2011; Foullon et al. 2011), the 
presence of magnetic field is not conducive to this instability. 
For example, magnetic field along the flow suppresses KHI 
whereas transverse (to the flow) field has no effect on the 
instability (Chandrasekhar 1961). However, magnetic field 
not only quenches shear instabilities but can also facilitate 
it. For example, the most important instability in accre- 
tion discs, the magnetorotaional instability is caused by an 
interplay between the Keplerian shear flow and magnetic 
field (Balbus & Hawley 1998). Thus, depending on the pres- 
ence of flow gradient, magnetic field, when well coupled to 
the plasma can suppress as well as drive shear instabilities. 
Therefore, before dwelling upon the role of magnetic field 
on the shear driven instabilities, it is pertinent to know how 
well is magnetic field coupled to the surrounding matter. 

The field drift through weakly ionised matter in the 
presence of shear flow can assist waves to grow. For example, 
in protoplanetary discs, both Hall and ambipolar diffusion 
assists magnetorotaional instability (Wardle 1999; Balbus 
& Terquem 2001; Kunz & Balbus 2004; Desch 2004; Wardle 
& Salmeron 2011; Pandey & Wardle 2012b). Clearly, drift 
of the magnetic field in a weakly ionised diffusive medium 
provides new pathways through which shear energy can be 
channelled to the waves. However, the role of diffusion in 
destabilising partially ionised medium in the presence of 
shear flow is not unique to weakly ionised discs but is quite 
generic (Kunz 2008; Pandey & Wardle 2012a). The crucial 
ingredients required to excite this diffusion-shear instabil- 
ity is: (a) the presence of a shear flow, and, (b) favourable 
magnetic field topology. The ensuing instability is overtly 
similar to KHI, and not surprisingly, the growth rate is pro- 
portional to the shear gradient. However, unlike KHI which 
is hydrodynamic in nature, this is magnetohydrodynamic 
instability. 

A detailed investigation of diffusive shear instability in 
the context of solar atmosphere is carried out in the present 
work which although similar to our previous work [Pandey 
& Wardle (2012a), hereafter PW12a] has following notable 
differences, (i) Unlike PW12a, where only vertical field and 
transverse fluctuation (vertical wavevector) is assumed, in 
the present case field topology is more general and wavevec- 
tor is oblique; (ii) In PW12a back reaction of the fluid on 
the magnetic field was completely ignored. By retaining the 
back reaction of the fluid on the field we show that shear 
driven diffusive instability do not have any cut-off wave- 
length; (iii) For vertical field and transverse fluctuations we 
showed in PW12a that only Hall diffusion assists the insta- 
bility and ambipolar and Ohm, which combines as Pedersen 
diffusion only causes damping of waves. In the present work 
for general field topology and oblique wavevector, we shall 
see that both ambipolar and Hall diffusion can assist the 
instability; (iv) The general stability criterion for magnetic 
diffusion dominated plasma in the presence of shear flows is 
presented in this work. Therefore, this paper is more general 
in scope and application than PW12a. 

The paper is organised in the following fashion. The ba- 
sic set of equations and dispersion relation is given in Sec. 2. 
In subsection 2.1 we give linearised equations in terms of dif- 
fusion tensor before describing general dispersion relation. In 
Sec. 3 the general stability criterion is described and maxi- 



mum growth rate of the instability is derived. The expression 
for maximum growth rate and critical wavelength is given 
in various limiting cases. In Sec. 4, various limiting cases of 
the dispersion relation is discussed. In Sec. 5 application of 
the result is discussed. Finally, in Sec. 6 brief summary of 
the result is given. 



2 BASIC MODEL 

The photosphere-chromosphere plasma consists of electrons, 
protons, singly ionized metallic ions, H, He I, He II, and He 
III. We shall ignore the distinction between hydrogen and 
metallic ions and assume that photosphere-chromosphere is 
made up of electrons, singly charged ions and neutral Hy- 
drogen. Although basic set of equations for partially ionised 
plasma was formulated more than 50 years ago (Cowling 
1957; Braginskii 1965), we shall use the single fluid formula- 
tion given by Pandey & Wardle (2008) for which the spatial 
and temporal scales of partially ionised plasma has been 
clearly elucidated. 

In the low frequency limit, the collisional dynamics and 
fractional ionisation can be incorporated in the basic set 
of equations without having to deal with the complexity of 
three-fluid equations. Thus the continuity equation of the 
bulk plasma fluid is 



! + v.(p„) = o. 



(1) 



Here p = pi + p n is the bulk mass density and pi : „ = 
nii ,„ m , n is the ion and neutral mass densities with 
the ion and neutral mass and number den- 
sities respectively; v — (pi Vi + p„ v„)/p is the bulk velocity, 
and, Vi and v n are bulk velocities of the ion and neutral 
fluids respectively. The momentum equation is 



at c 



(2) 



where J — en e (vi — v e ) is the current density, B is the 
magnetic field and P — P e + Pi + P n is the total pressure. 
The induction equation is 



~ {vxB) _ilL110j_^^j xb 

c c 

(j x b) x b 

c 



(3) 



where b = B/B is the unit vector along the magnetic field, 
and Ohm (rjo), ambipolar (tja) and Hall (tjh) diffusivities 
are 

„2 



D B 

no = -. — , va = - 

4 7TCT 4 7T Pi Vi, 

Here 



,VH = 



cB 
4 77 en e 



eerie 



Vi 



(4) 



(5) 



is the parallel conductivity, cj cj = eB/nij c is the particle's 
cyclotron frequency where e , B , rrij , c denotes the charge, 
magnetic field, mass and speed of light respectively and D = 
pn/ p is the ratio of neutral and bulk densities. For electrons 
v e = Ve n and for ions z/» = v in . Although v ee , v ei , vu and , v ie 
can become comparable to v en (see Table. 1), it is the 
neutral-plasma collision that gives rise to ambipolar and 
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Hall diffusion in the medium. 2 The electron-ion collision con- 
tributes to Ohm diffusion. 

Defining the plasma Hall parameter /3j as 

w, 



and the Hall frequency 



Pi 

P 



(6) 



(l) 



the diffsivities in Eq. (4) can be written in the compact form 
(Pandey & Wardle 2008) 



rjH 



,VA = D 



, and ,rjo = Pe 1 Vh , 



(«) 



where v ni = pt v in /p n and v A = B/^fTWp. 

The induction Eq. (3) can be explicitly written in terms 
of fluid and field velocities as (Wardle & Salmeron 2011) 



dB 

dt 



Vx 



{v + Vb)xB _^H°j^ 



where the field drift velocity vb is 



vb — r\p ■ 
and 



(V X B) x X b 
B 



(V X B) 



B 



rjp = rjo + rjA , 



(9) 



(10) 



(11) 



is the Pedersen diffusivity. The parallel and perpendicular 
components of current in the preceding equation refer to the 
orientation with respect to the background magnetic field. 
In the fluid frame, the magnetic drift velocity vb contains 
effect of field diffusion. As we shall see, the induction Eq. (9) 
in terms of magnetic drift velocity is particularly convenient 
for linearization which we will do shortly. 

We note from Eq. (4) that Ohm diffusion is independent 
of the magnetic field. Whereas Hall has linear dependence 
on the field, i.e. r\u oc B, the dependence of ambipolar diffu- 
sion on the magnetic field is quadratic, i.e. t\a oc B 2 . Thus, 
as the magnetic field strength in a flux tube decreases with 
increasing altitude, the drop in the Hall diffusivity will not 
be as severe as in ambipolar diffusion. The altitude depen- 
dence of diffusivities can be easily quantified if variation of 
the magnetic field with altitude is known. Thus, in order 
to compute diffusivities at various altitude, we shall choose 
power law variation of the magnetic field with neutral num- 
ber density n n 



B = B [ — 

. n 



(12) 



where no is the number density of neutrals at the surface 
(h — 0) and Bq = 1.2 kG is the typical value of the field at 
intergranular boundaries. Such a field profile captures essen- 
tial height dependence of observed field in flux tubes (Mar- 
tinez et al. 1997). We note that above scaling of the magnetic 
field differs somewhat from the mass density scaling (Leake 
et al. 2005) although we have retained the same power law 
index 0.3. 

The value of neutral number density and fractional ion- 
isation (Model C, VAL81), the ratio of neutral and bulk 

2 To leading order collisions between like plasma particles u ee , uu 
do not cause diffusion (Longmire & Roscnbluth 1956). 
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Figure 1. The Ohm (J70), Hall (t]h) and ambipolar (t]a) diffusion 
profiles in the photosphere-chromosphere are shown for (a) B = 
1.2 kG and (b) B = 120 G. 



densities D, and, various collision frequencies are given in 
Table 1. The collision cross-section is not known for the 
metallic ions-neutral collision, so we have employed the 
Messy-Mohr analytic formula for the cross-section (see ap- 
pendix). The resulting plasma- neutral collision frequency in 
the present work is an order of magnitude smaller than 
previously used values by Pandey et al. (2008) where a 
larger plasma-neutral collision cross-section was assumed for 
ion — H2 collision. 

We shall utilise the above magnetic field profile, along 
with the collision frequencies given in Table 1 to compute 
diffusivities which are given in Table 2. In Fig 1 diffusivi- 
ties are plotted against height which shows that Hall diffu- 
sion dominates Ohm and ambipolar diffusion in the photo- 
sphere and lower chromosphere in the intense field regions 
[Fig 1(a)]. However in the weak field (~ 100 G at h = 0) 
regions Ohm and Hall are comparable below 0.3 Mm before 
Hall becomes dominant diffusion in the entire photosphere- 
chromosphere [Fig 1(b)]. 

It has been suggested in the past that Pedersen diffu- 
sion (which is the sum of Ohmic and ambipolar diffusions) 
could play an important role in the chromospheric heating 
(Goodman 2004). We find that in the sunspots, active re- 
gions, pores and intergranular regions where field is intense, 
Hall and ambipolar will dominate Ohm diffusion. Clearly re- 
alistic modelling of various energetic processes such as chro- 
mospheric heating, CMEs etc. must include these diffusive 
processes. It would appear that since Hall diffusion is non- 
dissipative in nature, it may not have any role in the energy 
extraction from the convective motion of the largely neutral 
medium. However, as we have shown recently (Pandey & 
Wardle 2012a), Hall diffusion in the presence of shear flow 
destabilises low frequency fluctuations which in turn may 
lead to the turbulent cascade of convective energy to smaller 
scales where dissipation can convert it to heat. Therefore, it 
is quite plausible that the Hall effect may play an important 
role in heating the chromospheric plasma. 

How much diffusion is too much in the solar atmo- 
sphere? For example, if magnetic diffusion dominates fluid 
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Table 1. The neutral number density n n , the fractional ionisation X e = n e /n„ [Model C, VAL81], the 
ratio of neutral to bulk mass densities D = pn/p, the ion- neutral, electron-neutral, and electron— ion collision 
frequencies are given in the table. We have assumed B = Bo (n n /no)°' S with Bo = 1.2 kG and irii = 30 m p 
and m n = 2.3 m p where m p = 1.67 X 10 — 24 g is the proton mass. 
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Table 2. The values of Ohm, rjo, ambipolar, r\A and Hall, rjrj diffusivities at different altitude are given in the table. The value of 
various collision frequencies have been taken from Table 1. 
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convection in the induction Eq. (3), field will be poorly cou- 
pled to the plasma (Wardle 2007). Thus, we shall compare 
advection term Vx(t)XfJ)~«B/L with various diffusion 
term in the induction Eq. (3) by defining magnetic Reynolds 
numbers Rm , Am and Hm as 

Rm = ^ )Am= V JL &Hm= — . (13) 
rjo Va Vh 

Note that both Am and Hm depends on how well the plasma 
is coupled to the magnetic field since both ambipolar and 
Hall diffusion depends on the magnetic field. Further tja = 
Pi t)h and r) H = /9 e J?o • Thus 

Am = Rm/ (ft /3 e ) , Hm = Rm/ft . (14) 

The dependence of Am and Hm on the plasma Hall param- 
eter /3j is not surprising since the origin of both ambipolar 
and Hall diffusion is inherently linked to the magnetisation 
of the medium. This also explains inherently different nature 
of the Ohm and ambipolar diffusion: whereas Ohm diffusion 
acts isotropically, ambipolar diffusion owing to its depen- 
dence on the magnetic field is anisotropic. As we shall see, 
in the presence of shear flow, the anisotropic nature of am- 
bipolar diffusion is at the centre of wave destabilization. 

The height dependence of Rm, Hm and Am for both kG 
and 0.1 kG fields were discussed in PW12a. It was shown 
that when v ~ va (where va is the Alfven speed), Rm 3> 1 
suggesting that the Ohm diffusion is unimportant in com- 
parison with the advection term. However, Hm and Am are 
three to four orders of magnitude smaller than Rm. The 
recent 2D numerical simulation of partially ionized solar at- 
mosphere also suggests that in the weak field (< 100 G) 
regions in the chromosphere, the Hall diffusion is two or- 
ders of magnitude larger than the Ohm diffusion whereas 
ambipolar diffusion is four to six order of magnitude larger 



than the Ohm diffusion (Sykora et al. 2012). Clearly, both 
ambipolar and Hall diffusion terms are important in the in- 
duction equation. 



2.1 Dispersion relation 

The magnetic element in the photosphere-chromosphere 
which can be modelled as cylindrical tubes or planer sheets 
are highly dynamic accompanied by numerous flows with 
different spatial/ temporal scales. Recent numerical simula- 
tions of the umbral magneto convection suggests that the 
dynamical scale over which Hall diffusion can generate mag- 
netic and velocity fields is much faster [~ 10-20 km spatial 
scale and ~ 300 s temporal scale (Cheung & Cameron 2012)] 
than the spatial and temporal scale of a typical flux tube (> 
few hundred km and ~ few days). The simulation results 
are easily scalable to 2 km with temporal scale ~ 2 s. We 
note that at present the best achievable resolution is 90 km 
(Bonet et al. 2008). 

Since the spatial scale over which the flow and field 
generation occurs, is much smaller than the typical tube 
diameter, we shall approximate the cylindrical tube by a 
planer sheet and work in the Cartesian coordinate system 
where x , y , z represents radial, azimuthal and vertical di- 
rections locally. We shall assume an initial homogeneous 
state with azimuthal shear flow v — vo' xy. The magnetic 
field in the intergranular lanes at the network boundaries is 
clumped into elements or flux tubes that are generally ver- 
tical (Martinez et al. 1997; Hasan 2009) but highly inclined 
fields have also been reported in the literature (Stenflo et 
al. 1987; Solanki at al. 1987). The internetwork magnetic 
elements have predominantly horizontal field (Hasan 2009; 
Steiner & Rezaei 2012). Therefore, we shall assume uniform 
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background field that have both azimuthal as well as vertical 
component, i.e. B — (0, B y , B z ). 

The focus of the present investigation is low frequency 
behaviour of the medium, and thus, we shall work in the 
Boussinesq approximation limit which is valid if the mo- 
tion in the medium is very slow (Spiegel & Veronis 1960). 
Thus linearising Eqs. (1) and (2) and assuming an ax- 
isymmetric perturbations of the form exp (ik ■ x + at), with 
fc = (k x ,0, k z ), in Boussinesq approximation, we get 



k ■ Sv = . 



a dv x = — i k x h 



(15) 



p 4n p 

a Sv y + v' Sv. 
8p i 



(k ■ B) SB X - (B ■ SB) k x ] , 

(fc ■ B) SBy , 



a5v z = -ik z — + 

p 4 7T p 



4n p 

l(k-B) 8B z -(B-SB)k z ] . (16) 



Eliminating the pressure perturbation in favour of ve- 
locity and making use of Eq. (15), from preceding equation 
we get for the (x , y) components 



a 2 \ -v a J 



(17) 



Here Sv = Sv/va and SB — SB/B , v' = dv(x)/dx, p = 
^fc ■ b^j = k z b z , fc = k/k and k z — k z /k. Since 

Sv B =ik ^rjH |/i SB x b - (b ■ SB^ fc x b j 

+ VP ^SB-(b-SBj fc}] , (18) 

the linearised induction equation (x , y components) can be 
written as 



G 

-v'o 



l 2 2 2 
k V A p 



a 

-v'q a 



+k 2 



I T/ Xx 


Tjxy \ 


SB 


V r lv x 


Vyy ) _ 





(19) 



sm +9VA 

') 



Tjxy — 

no + (l-k 2 x b 2 z ) n, (20) 



where rj is the diffusivity tensor with following components 

r/xx = r/o + b 2 z r/A 
Vyx = (gVA- st) H ) /k 2 z , r) yy -- 
Here 

g = -kx k z b y b z , (21) 

and helicity s — p k z = b z k 2 . We note that for purely verti- 
cal field the sign of helicity s is determined by the projection 
of the vertical magnetic field on the vorticity, Vxi>o- 

Following dispersion relation can be derived from (19) 



a 4 + k 2 (t7p + rj T ) cr 3 + C 2 o 2 + Ci a + Co = , 
where 

C 2 = k 4 (rj P rj T + p 2 r)li) +k 2 [2 p 2 - a (g fj\ + s r/ H )] 



(22) 



a = -v'q/ \v' \ 
Vt = no + /J 2 r\ A . 



Ci = k p (-Pip + r) T ) , 
Co = fc 4 p 2 [p 2 - a (gr) A + sfju)] , (23) 
±1 and 



(24) 



We have used following normalisation in the above equations 
a-— k- kVA f,-" 1 ^ 

bo I Kl 



(25) 



Eq. (22) reduces to known dispersion relation (Kunz 2008) 
when D = 1 and rjo — 0. When D = 0, ambipolar diffusion 
expectedly, drops out of the dispersion relation. 



3 NON-IDEAL MHD INSTABILITIES 

A general stability criterion in the magnetic diffusion dom- 
inated plasma can be derived from the dispersion relation, 
Eq. (22) by recasting it in the following simple form 



a fc 4 + b fc 2 + c : 



0, 



where the coefficients a, b and c are 



jp 



b = 8- 



~2 

■ya . 



(26) 
(27) 



a = P 
with 

P = {tJ 2 vh + vp ffr) 5- 2 + (vp + m) P 2 o + p 4 , 

8 = [{ftp + r) T ) o + 2 p 2 ] a 2 , 

7 = a (gfjA + s?7h) • (28) 
The dispersion relation, Eq. (26) can also be recast as 

Pk 4 + Sk 2 + & 4 = 7 fc 2 \a 2 + p 2 fc 2 ] . (29) 

Thus in the neighbourhood of a ~ fc <JC 1 where 

p ~ p 4 , 8 ~ 2 p 2 a 2 , (30) 

Eq. (29) becomes 



2 72 . -2 

p k + a 



2 ~2 

7 a , 



which can be written as 

2 



7-/i 



(31) 



(32) 



From Eq. (32) we see that for positive a right hand side is 
positive. This implies that 7 — p 2 > 0. Therefore, we arrive 
at the general stability criterion which states that if 



a (gVA + sf/H) > p , 



(33) 



the waves are unstable in the medium. The Ohm diffusion 
do not appear in the above expression which is not sur- 
prising considering that the above criterion pertains to the 
long wavelength fluctuations. For definiteness, in subsequent 
analysis, we shall assume a = 1. 

The diffusion-shear instability is caused by a competi- 
tion between the fluid advection and field drift in the plasma. 
This can be seen from the y and x components of Eqs. (17) 
and (18) respectively, which suggest that when a = 0, the 
advection of fluid in the x direction, 8v x is equal and oppo- 
site to the magnetic field drift velocity Svb x - By combining 
y component of Eqs. (17) and x component of Eq. (18)] we 
get 

.2 



Sv x + Svbx =ikp 



^f- SBy + l- fj A SB, 

L ap 2 k 2 



(34) 



Here we have neglected Ohm diffusion. By setting — ~f + p 2 = 
near marginal state, in the vicinity of tja <S 1, we get 
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5v x + Svbx ~ 0. This provides a simple physical explana- 
tion on how the magnetic diffusion helps the shear flow in 
destabilising the waves. The outward slippage of the field 
in the x-direction due to magnetic diffusion weakens the 
magnetic tension force. As a result magnetic restoring force 
in the wave is dominated by the inertia force, resulting in 
the increased inward drift of the fluid elements. This is how 
magnetic diffusion assist the waves to grow in the presence 
of shear flow. 

The stability criterion, Eq. (33) in the dimensional form 
becomes 



-v' (gr/A + stjh) > v 2 



(35) 



which suggests that when — v' > 0, the linear combination of 
ambipolar and Hall diffusivities multiplied by suitable topo- 
logical factors g and s must exceed square of the oblique 
Alfven speed. Above equation provides a simple stability 
criterion of diffusive medium. For example, low frequency 
fluctuations in the Hall-Ohm dominated photosphere-lower 
chromosphere (r/A = 0) are unstable if 



-v'qVh > b z v 2 A , 



(36) 



or, in terms of Hall frequencies 

—v'o > b z ujh = uhz ■ (37) 

Here we have used t)h = v\/u>h (PW12a) and b z ujh = ujhz 
is the Hall frequency defined in term of vertical field. Clearly 
if the Hall frequency is less than the shear frequency — v'o, 
Hall diffusion can drive shear flow instability. 

When v'oTih/va = —b z , since t\a > 0, ambipolar dif- 
fusion can drive shear flow instability provided g > 0. 
For purely vertical field or, transverse fluctuations (vertical 
wavevector) when g = 0, ambipolar diffusion can only cause 
the damping of waves. Therefore, the question of ambipo- 
lar diffusion assisting the shear flow instability is inherently 
linked to the ambient field geometry and wave obliqueness 
which together is encapsulated in the topological factor g. 
The important role of g in the ambipolar diffusion driven 
shear flow instability was discovered by Desch (2004). 

How does ambipolar diffusion help drive shear flow in- 
stability? In order to see this, we first note that when the 
wavevector is perfectly aligned to the ambient magnetic 
field, i.e. fj, = k z b z = 1, ryr = rjp in Eqs. (22)-(24), Ohm and 
ambipolar diffusion combines together as Pedersen diffusion 
and their effect on the wave propagation is identical-they 
both cause wave damping. Only when n / 1, this diffusive 
degeneracy is lifted and Ohm and ambipolar diffusion can 
no longer be combined together as single diffusion. After 
removal of degeneracy whereas Ohm diffusion still causes 
isotropic damping of the waves, damping by the ambipolar 
diffusion becomes anisotropic. Therefore, the non-vertical 
magnetic field and oblique wavevector (which is crucial in 
removing this degeneracy) plays important role in the am- 
bipolar diffusion driven shear instability. This can also be 
seen from the genelarlised induction Eq. (3) if we rewrite 
the electric field as 

E' = ^ { Vo + [1 - (j ■ b) 2 ] VA } J + ^-VhJ X b 

+ r !A (j ■ b) [j X (J X b)] , (38) 

where the electric field is written in the neutral frame, j — 
J/\J\- It is clear from the preceding equation that when 



medium is threaded only by the vertical field and wave is 
propagating along the field, i.e. k — (0, 0, k z ), j = (jx,jy,0), 
both Ohm and ambipolar diffusion cause damping of the 
waves since (j ■ b) = and E' ■ J — rjp J 2 . However, when 
/j, 7^ 1, last term in Eq. (38) can help fluctuations grow since 

E ' ■ J = ^ {no + [i - U ■ b) 2 ] va } J 2 



+ ^-t,a (j-b) [J ■ (j X (J X 6))] 



(39) 



Therefore, ambipolar diffusion plays dual role in a partially 
ionized medium: whereas, in one direction it can cause dis- 
sipation like Ohm, in the other direction the dissipation is 
considerably smaller. The directional dissipation is the hall 
mark of ambipolar diffusion which causes it to assist the 
instability (Desch 2004). 

From stability criterion Eq. (35) it is clear that the 
Hall-Ohm stable medium (— v' 7]h < b z v\) may or may not 
be ambipolar unstable whereas Hall-Ohm unstable medium 
(—v'o rjH > b z Va) can always become ambipolar unstable for 
non-vertical field and oblique wavevector since g > can 
be easily satisfied. The ambipolar diffusion not only drives 
shear flow instability when g > 0, but also enlarges the 
parameter space over which Hall can destabilise the waves. 
Therefore Hall diffusion can drive shear flow instability when 



(40) 



Since b z ,b y G [0, 1] and maximum g — 0.25, above inequal- 
ity implies that for non-zero tja, Hall diffusion drives shear 
instability at much larger negative value than when the field 
has only vertical component. 

The stability criterion can also be recast as 



kx by 
T z b z 



> 



r/A 



b z 



(41) 



Assuming positive left hand side in the preceding equation 
we see that above criterion is easily satisfied for non-zero r/A 
when k z — > 0. Therefore, waves propagating almost along x 
direction i.e. when the fluctuations are almost magnetosonic 
are always unstable. 

The maximum growth rate of the instability can be 
found by setting the discriminant b 2 — 4 a c = in Eq. (26). 
This yield 



gr]A + sr/H 



[(??p + tit) + 2 ^/n 2 r 1 2 H + qp rjrj 



(42) 



We see that the maximum growth rate depends on both g r/A 
and s t)h and for comparable r/A and rjH [as is the case in 
large part of the solar atmosphere; see Fig. (1)], the bigger 
contribution to the growth rate comes from the Hall diffusion 
since the maximum value of parameter g is 0.25 whereas 
maximum value of the helicity s is one. 

In the absence of Hall and Ohm, the maximum growth 
rate, Eq. (42) becomes 



(Jo 



9 Pol 



(43) 



(i + m) 2 ' 

which suggests that when the field drift is solely due to the 
ambipolar diffusion the maximum growth rate is indepen- 
dent of the diffusivity as well as the strength of the back- 
ground magnetic field. However, the signature of the mag- 
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(a) B = 1.2kG 
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Figure 3. Above cartoon depicts growth rate against wavcnum- 
ber. 



Figure 2. The maximum growth rate, Eq. (42) is shown for k z = 
b z = 1 (bold line), k z = k x = b z = — b y = 1/V2 (dotted line) and 
k z = 1/^2 and b z = 0.1 (dashed line) for (a) B = 1.2kG and (b) 
B = 120 G. 



netic field in Eq. (43) appears through the topological factors 
g and /i. 

When the field drift in the plasma is solely due to the 
Hall diffusion, the maximum growth rate Eq. (42) becomes 



g ??A + s fju > 2 + 



fjp + J?t 



<T = 



2 



(44) 



which implies that the waves with k z = 1 are most unstable. 
A comparison with Eq. (43) shows that in purely ambipo- 
lar or, purely Hall case, growth rate is independent of the 
ambient diffusivity. 

We see from Fig. 2(a) that for purely vertical field and 
vertical wavevector (k z = b x = 1, g = 0, bold lines) the in- 
stability grows at maximum rate |i>o|/2 in the photosphere- 
lower chromosphere in the strong field region. Recall that 
in this interval Hall is the dominant diffusion mechanism in 
the network and internetwork regions (Fig. 1) and dissipa- 
tion due to the Ohm and ambipolar diffusion (which can be 
combined together as Pedersen diffusion for this topology) 
is small. 

When the field is weak [Fig. 2(a)] the instability 
grows at close to maximum rate in the entire photosphere- 
chromosphere except very close to the surface (< 0.2 Mm). 
When k z — k z — b z — ~b y — l/\/~2 [dotted curves in 
Fig. 2(a)-(b)] the growth rate is smaller than the previ- 
ous case. This is because in this case the instability is not 
only due to the Hall diffusion but is also assisted by the am- 
bipolar diffusion through directional dissipation of waves. 
Therefore, the growth rate is always smaller when both Hall 
and ambipolar diffusion are present in the medium. With 
the decreasing vertical field [b z — 0.1, dash-dot curves in 
Fig. 2(a)-(b)] the maximum growth rate decreases further 
owing to smaller g and s. 

The most unstable wavenumber is 



(45) 



M Vh +fpfr cr — fi 



It is clear from preceding equation that when 



(46) 



VP rjT 

maximum growth rate occurs at finite k [see Fig. (3)]. 
In the opposite limit fco becomes imaginary implying that 
maximum growth rate occurs at infinity. For parameters of 
network-internetwork magnetic elements, with shear gradi- 
ent v'q ~ 10~ 2 s _1 and va ~ 10 5 kms -1 , fco is imaginary 
[corresponding to two similar curves in Fig. (3)]. 



4 LIMITING CASES OF THE DISPERSION 
RELATION 

The magnetic field drift in the partially ionized plasma opens 
up new pathways through which the free shear energy of the 
fluid can be transferred to the waves (PW12a). To see this, 
we first note that since a ~ 1, three terms in Eq. (19) are 
~ 1 , k 2 , k 2 rj and thus, one of the following three scenarios 
may prevail in a diffusive medium. 

A. Ideal MHD: In this limit k 2 ~ 1 > k 2 r), i.e. f] < 1 
and last term in Eq. (19) can be neglected. However, as 
has been shown in PW12a, this limit is applicable to the 
long wavelength fluctuations > 10 3 km and therefore, is not 
relevant to the present analysis. 

B. Cyclotron limit: In this case k 2 f] ~ k 2 3> 1. This 
is the low frequency limit and first term in Eq. (19) can 
be neglected. The low frequency, short wavelength dressed 
ion-cyclotron wave with frequency uic = M ^>h is the normal 
mode of the system (Pandey & Wardle 2006, 2008). 

Assuming \v' \ ~ 0.01s -1 and va ~ 5 x 10 5 kms" 1 we 
get fj = v\ v 'o\/ v a ~ 10~ 2 — 10~ 4 <g; 1 and thus it was in- 
ferred (PW12a) that the cyclotron limit is not valid in the 
photosphere-chromosphere. However, at an increased shear 
frequency \v' Q \ < Is -1 we get fj = ri\v' Q \/v 2 A < 1 and this 
limit becomes important. Since recent numerical simulation 
with 10 — 20 km resolution (Cheung & Cameron 2012) can 
be easily scaled to 2 km which for typical vo ~ 2 km s - gives 
| Wo | ~ Is 1 ) we conclude that the proper analysis of the cy- 
clotron limit will provide important insight to the ongoing 
numerical simulations of the photosphere-chromosphere. 

Taking k — > oo, and thus setting a(er) = in Eq. (26) 
or, neglecting first matrix in Eq. (19), we get following dis- 
persion relation 

(m 2 m. + Vp Vt) 



a 2 + ji 2 {fjp + f/ T ) cr = jj, 2 (7 — jj, 2 



(47) 
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Positive a requires 7 > fi 2 since coefficients on the left hand 
side of preceding equation are positive. Thus the stability 
criterion in the cyclotron limit is identical to general sta- 
bility criterion, Eq. (33). This is not surprising since above 
dispersion relation is short wavelength limit of the general 
case. 

The growth rate of the instability becomes 



2 (ijp ffr + M 2 Vh) 



[- (fjp + ffr) + Va] 



where 



A = (fj P - ffr) 2 + 4 rjn (7 - /j 2 ) + 4 



7^p Vt 



(48) 



(49) 



Above equation acquire particularly simple form in purely 
Hall or ambipolar limits. For example in the Hall diffusion 
dominated regime, from Eq. (48) we get 



a = 



1 s 2 va 

-v /X — 

1JH Vh 



1/2 



VA : 



which in ix < 1 limit becomes 

~ r ' 1 1/2 
a w [-v soj h \ 



(50) 



(51) 



Thus the growth rate of the dressed ion-cyclotron wave ap- 
proximately equals the geometric mean of the shear and Hall 
frequencies and attains maximum value only for positive he- 
licity s = 1. 

In purely ambipolar diffusion dominated case setting 
no =t)h = 0, we get 



A. 

2va 



^gv'oVA 



v\ 



(52) 



We see from Eq. (52) that when g — 0, ambipolar diffusion 
causes only damping. As noted in the previous section, it is 
only when topological factor g is non-zero that anisotropic 
ambipolar diffusion can drive shear flow instability. 



C. Highly diffusive limit : In this limit k fj 



1 > 



fc 2 , i.e. fj ^> 1 and fc 2 <JC 1. Only the first and last 
term in Eq. (19) are retained. For typical values of 7/ [Ta- 
ble (2)], this limit gives A < 6 km which fits within a 
pressure scale height. Therefore, as noted in our previous 
work (PW12a), highly diffusive limit is applicable to the 
photosphere- chromosphere. 

The dispersion relation in this limit becomes 

(fjp + fj T ) ak 2 + a 2 = k 2 [7 - fc 2 (fj P fj T + /J 2 ??h)] , (53) 

from where we see that for positive a left hand side is posi- 
tive. Thus right hand side must be positive. Thus we arrive 
at following general stability criterion 

-v'o (gVA + sr/ H ) > k 2 (fjp fj T + fj 2 fin) . (54) 
In the absence of Hall diffusion, above criterion becomes 



-v' grj A > k 2 rj P rj T ■ 



(55) 



As right hand side in the above equation is positive, the 
stability criterion implies that in the presence of favourable 
shear flow gradient (— v' > 0) ambipolar diffusion will drive 
the instability only if g > 0. Positive topological factor g 
guarantees that the ambipolar diffusion will assist the shear 
flow in destabilising the waves. 



When ambipolar diffusion is negligible, Hall diffusion 
can as well drive shear instability provided 

72 l~2 , 2 -2> 



-v'o srjH > k 2 (r?o + m) 



(56) 



Above criterion is similar to Eq. (24) of PW12a. Note that 
in PW12a, s = 1 as both the field and the wavevector is 
vertical whereas here we are dealing with the general field 
topology (s 1). 

The instability in the Hall-ambipolar diffusion domi- 
nated regime will grow at the maximum rate 



(J = 7 



(jjp + fjr) - 2 ^/ffp ifr + [i 2 ff^ 
(fjp -vt) 2 -4/j 2 r]|, 



(57) 



which for purely Hall (t)a — r\o — 0) case reduces to 
Eq. (44). In the absence of Hall diffusion the maximum 
growth rate becomes 



9va 



(58) 



When rjo = above equation reduces to Eq. (43) . 

Most unstable wavelength in the highly diffusive limit 
becomes 



k 



2a 2 



[a (grjA + s?7h) - (r~/p + t]t) &o] 



(59) 



For purely Hall case and (a — 1) preceding acquires partic- 
ularly simple form 



fco 



1 



2 b z fj-R 

In purely ambipolar regime fco becomes, 
1 



fco = 



(1 + M) V M?7A ' 



(60) 



(61) 



from where it is clear that the most unstable wavenumber is 
nonzero only when the topological parameter g is nonzero. 
This is not surprising given that very existence of the am- 
bipolar diffusion driven shear instability depends on the field 
geometry and obliqueness of the wave vector. 



5 DISCUSSION 

The solar photosphere is threaded by kG field concentrated 
in the vertical flux tubes (radius ~ 100 — 200 km) at in- 
tergranular boundaries (Hasan 2009). Although less fre- 
quently (< 40%), field of similar strength have also been 
observed in the quiet solar internetwork region (Almeida et 
al. 2010). Outside kG patches, internetwork field is much 
weaker (~ 100 G). The flux tube is often modelled as non- 
rotating cylindrical tubes with the plasma and the magnetic 
pressure balance providing required stability. In the present 
work we have approximated flux tubes by a planer sheet 
where x and y coordinate locally correspond to the radial 
and azimuthal directions locally. 

Although the chromosphere has limited extension 
in comparison to the corona, its net radiative loss ~ 
10 7 erg cm -3 s _1 is 10 times more than the corona. Further, 
except for the flares, most solar atmospheric heating occurs 
in the chromosphere (Aschwanden 2001; Goodman 2001). A 
strong correlation between core emission of calcium K and 
H resonance lines and quite sun magnetic field (Schrijver et 
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al. 1989) suggest that the origin of chromospheric heating 
(~ 10 7 ergs cm" 2 s~ x ) is magnetic in nature. The magnetic 
diffusion driven shear instability proposed in the present 
work can provide a viable mechanism to explain the excess 
chromospheric heating since crucial ingredients required to 
excite this instability namely shear flow and magnetic field 
is always present in the network-internetwork region. The 
attractive feature of this fast growing diffusive shear instabil- 
ity is that all wavelengths of fluctuations are likely to be ex- 
cited as there is no cut-off wavelength. Thus, as we see from 
Fig. (2) , network field below 1 Mm is likely to be destabilised 
by this instability. In the internetwork elements, this insta- 
bility may operate in the entire photosphere-chromosphere. 

Only uncertainty involved is the lack of information 
about the scale of shear flow gradient which can not be re- 
solved by current observations. However, small whirlpools 
with size similar to terrestrial hurricanes [< 0.5 Mm with 
typical lifetime ~ 5 min. on the solar surface Bonet et al. 
(2008)] suggest the presence of such flows. Long lasting large 
scale vortices at supergranular junctions with typical life- 
time ~ 1 — 2 h with enhanced Call emission have also been 
observed (Attie et al. 2009). The swirl motion in the chro- 
mosphere has been recently detected by Wedemeyer-Bohm 
& Voort (2009). Clearly, observations suggest ubiquitous 
presence of flow gradients in the photospheric-chromospheric 
plasma. Past numerical simulations also indicated the pres- 
ence of vortex flows in intergranular lanes (Zirker 1993; 
Stein & Nordlund 1998). The typical vorticity of a vor- 
tex is ~ 6 x 10 -3 s _1 which corresponds to rotation period 
~ 35 minutes (Bonet et al. 2010). Thus it would appear 
that Hall instability does not have time to develop since the 
growth rate (|uq|/2 = 3 x 10~ 3 s" 1 ) is very small. However, 
above vorticity value is limited by the upper limit in the 
vorticity resolution [~ 4 x 10~ 2 s -1 , Bonet et al. (2010)]. 
The numerical simulation gives much higher vorticity value 
(~ 0.1 - 0.2 s 
(Fig. 31, Stein 
sponding to \v' \ = 0.2 s" 1 is one minute. 

For almost magnetosonic waves (k z — > 0), maximum 
growth rate of the Hall diffusion driven shear instability may 
becomes quite small. The maximum growth rate in the am- 
bipolar diffusion dominated middle and upper chromosphere 
will be only one fifth of the Hall diffusion dominated case 
for the maximum g — 0.25 and n — 0.5. Therefore, vortex 
motions with typical lifetime > 15 min. will be susceptible 
to the ambipolar diffusion driven shear instability. Since, 
vortex motions of various spatial and temporal scales are 
observed, it is likely that non-ideal MHD effects will play 
important role in exciting low frequency turbulence in the 
medium. 

The non-ideal MHD description of the photosphere- 
chromosphere provides several pathways through which 
shear energy can be channelled to the waves by magnetic 
field. For example in excessively diffusive limit when Sv /v <C 
SB/B, diffusion in tandem with the shear flow can desta- 
bilise the network-internetwork field. For purely vertical field 
and vertical wavevector this limit has been discussed in de- 
tail in PW12a. In order to compare with PW12a, we briefly 
describe the effect of the field topology and wave orienta- 
tion in the highly diffusive limit. Comparing Figs. 2 and 4 
we conclude that the maximum growth rate is similar in 
both cases. When b z = k z — 1, instability grows at maxi- 



B = 1 .2 kG 
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; Nordlund (1998)). The growth rate corre- 




1 1.5 
Height [Mm] 

B = 120G 



0.5 



(a) 

- ■ • ■ mm*''/, 











0.5 



1 1.5 
Height [Mm] 



2.5 



0.25 
0.2 



«:° 0.1 





k /K = 1,B Z /B = 1 




■ ■ »k /k=0-9, B ; / B = 0.9 






mmmm\ik = QA, B ; / B = 0.9 









1 1.5 
Height [Mm] 



Figure 4. The growth rate and most unstable wavelength is 
shown for 1.2 kG [Fig. 4(a)] and 120 G [Fig. 4(b)] fields. Following 
parameters have been used in the above figure: k z = 1 , b z = 1 
(bold line), k z = b z = 0.9 (dashed line) and k z = 0.4, b z = 0.9 
(dotted line). 



mum rate whereas with decreasing k z or, b z the growth rate 
diminishes. 

In Fig. 4(b) we plot most unstable wavelength against 
height for same parameters as in Fig. 4(a). The wavelength is 
normalized against scale height calculated self-consistently 
using model C, VAL81. For both kG and weaker field Ao fits 
well within a scale height and thus the instability will grow 
at maximum rate in the entire photosphere-chromosphere. 
However, it is only in the photosphere and lower chromo- 
sphere (< 1 Mm) where the instability will grow at maxi- 
mum rate for a kG field. In the middle and upper chromo- 
sphere the growth rate tapers off and becomes one eighth 
of shear frequency. Therefore, in the strong field region, dif- 
fusive instability will be efficient in destabilizing magnetic 
element in the photosphere and lower chromosphere. When 
the field is weak (~ 100 G) , the instability can operate in 
the entire photosphere-chromosphere at maximum rate. 

How does nonlinear saturated state of diffusive insta- 
bility will look like? The answer to this can be given only 
by numerical simulations. However, if the nonlinear results 
of the protoplanetary discs and star forming regions are any 
guide then this instability should be quite efficient in excit- 
ing the low frequency turbulence and heating of the plasma. 
In fact interplay between the vortex flow and magnetic dif- 
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fusion could be responsible for entire energy budget of the 
solar atmosphere. 

The energetic events such as hard x-ray emissions (~ 
10 26 ergs -1 ) are believed to be due to the presence of ener- 
getic electrons with energies above (~ 20keV) in the solar 
flares. One of the plausible scenarios for the electron ac- 
celeration to such high energies requires turbulent medium. 
The stochastic acceleration of the electron by the turbulent 
wave spectrum produces high energy spectrum (Fletcher & 
Hudson 2008). Since magnetic diffusion could be important 
agent in driving the diffusive shear instability, this could 
easily lead to the low frequency turbulence in the medium. 
However, further work is needed in this direction to support 
this hypothesis. 



6 SUMMARY 

The granular motion is responsible for the generation of low 
frequency Alfvenwave in the predominantly neutral photo- 
sphere. Unlike high-frequency MHD counterpart which un- 
dergo damping (Vranjes et al. 2007, 2008), low frequency 
Alfvenwave generally propagate undamped in the medium. 
In the presence of shear flow gradient, the solar photosphere- 
chromosphere region can become unstable due to non-ideal 
MHD effects. Depending on the shear scale, such instability 
can be excited at all wavelengths. The presence of favourable 
magnetic field topology will facilitate the transfer of shear 
energy to the magnetic fluctuations. The instability in the 
solar atmosphere is due to the presence of two unrelated 
physical processes (a) nonideal MHD effects and, (b) shear 
flow. 

The following are the itemized summary of the present 
work. 

1. The magnetic diffusion in the presence of shear flow 
makes solar atmosphere susceptible to diffusive shear insta- 
bility. The instability depends on the sign of the shear gra- 
dient and on the local magnetic field topology. 

2. Purely vertical magnetic field and vertical wavevec- 
tor along with the shear flow is fundamental to this Hall 
diffusion driven instability. 

3. Only when field has non-vertical component and 
waves are propagating obliquely, both Hall and ambipolar 
diffusion together can assist this instability. 

4. The maximum growth rate of the instability which 
is proportional to the absolute value of the shear gradient 
occurs when both the field and the wavevector is vertical. 

5. The e-folding time of the magnetic diffusion driven 
shear instability is very short (~ 20 s for v' = 0.1 s _1 ) sug- 
gesting that the flux tubes are likely to be susceptible to this 
fast growing instability. 



ACKNOWLEDGMENTS 

The financial support of Australian Research Council 
through grants DP0881066 and DP12010I792 is gratefully 
acknowledged. This research has made use of NASAs Astro- 
physics Data System. 



REFERENCES 

Almeida J. S., Bonet, J. A., Viticchie, B. & Del Moro, D., 

2010, ApJ, 715, L26 
Almeida J. S., Gonzalez, M. M., 2011, as- 

troph/arXiv: 1 105.0387vl 
Arber T. D., Haynes, M. & Leake J. E., 2007, ApJ, 666, 

541 

Aschwanden, M., 2009, Phyics of the Solar Corona, (UK: 

Sringer-Praxis) 
Aschwanden, M. J. 2001, ApJ, 560, 1035 
Attie R., Innes D. E. & Potts, H. E., 2009, A&A, 493, L13 
Balbus S. A. & Hawley J. F., 1998, Review of Mod. Phys., 

70, 1 (BH98) 
Balbus S. A. & Terquem C, 2001, ApJ, 552, 235 
Balmaceda, L., Vargas Dommguez, S., Palacios, J., Ca- 

bello, L, & Domingo, V. 2010, A&A, 513, L6 
Bonet, J. A., Marquez, I., Almeida, J. S., Cabello, I. & 

Domingo, V., 2008, ApJ, 687, L131 
Bonet, J. A., Marquez, I., Almeida, J. S., Palacios, J., Pll- 

let, V. M., Solanki, S. K. et al., 2010, ApJ, 723, L139 
Braginskii, S. I. 1965, Review of Plasma Physics (vol 2) 

ed. Leontovich, M, A., 1965, 205 (New York: Consultants 

Bureau) 

Chandrasekhar, S. 1961, Hydrodynamic & Hydromagnctic 
Stability, (Cambridge: Cambridge University Press) 

Cheung, C. M. M. & Cameron R. H., arXiv: 1202. 1937 

Cowling, T. G. 1957, Magnetohydrodynamics, (Bristol: 
Adam Hilger) 

Desch, S. 2004, ApJ, 608, 509 

de Wijn A. C, Stenflo, J. O., Solanki, S. K. & Tsuneta S., 

2009, Space Sci. Rev., 144, 275 
Erdelyi R. & James S. P., 2004, A&A, 427, 1055 
Fletcher L. & Hudson H. S., 2004, ApJ, 675, 1645 
Fontana, P. R., 1961, Phys. Rev, 123, 1865 
Foullon C, Verwichte, E., Nakariakov, V. M., Nykri, K. & 

Farrugia C. J., 2011, ApJ, 729, L8 
Glassgold, A. E., Krstic, P. S.& Schultz, D. R., 2005, ApJ, 

621, 808 

Goedbloed, H. & Poedts S. 2004, Principles of Magnetohy- 
drodynamics with Applications to laboratory and Astro- 
physical Plasmas, (London: Cambridge) 

Goodman, M. L. 2001, Sp.Sci.Rev, 95, 79 

Goodman, M. L. 2004, A&A, 416, 1159 

Goodman, M. L. 2011, ApJ, 735, 45 

Hasan, S. 2009, in N. Goplaswamy & D. F. Webb (eds.) 

Universal Heliophysical Processes, IAU Symp. 257, 121 
Hasan, S. & Kalkofen, W. 1999, ApJ, 519, 899 
Khodachenko, M. L., Arber, T. D., Rucker, H. O. & 

Hanslmeier, A., 2004, A&A, 279, 389. 
Khomenko E. & Collados M., 2012, ApJ, 747, 87 
Karpen, J. T., Antiochos, S. K., Dahlburg, R. B. & Spicer 

D. S., 1993, ApJ, 403, 769 
Kolesnikov, F., Bunte, M., Schmitt, & Schussler M., 2004, 

A&A, 420, 737 
Krstic, P. S.& Schultz, D. R., 2009, Phys. Plasmas, 16, 

053503. 

Kudoh, T. and Shibata, K.: 1997, in A. Wilson (ed) Proc. 

of the Fifth SOHO Workshop, ESA SP-404, pp. 477 
Kunz, M. W. & Balbus S. A., 2004, MNRAS, 348, 355 

(KB04) 

Kunz, M. W., 2008, MNRAS, 385, 1494 (K08) 



© 2012 RAS, MNRAS 000, 1-?? 



12 B.P.Pandey and Mark Wardle 



Leake J. E. & Arber T. D., 2006, A&A, 450, 805 

Leake J. E., Arber T. D. & Khodachenko M. L., 2005, 

A&A, 442, 1091 
Lites B. W., Kubo M., Socas-Navarro H. et al., 2008, ApJ, 

672, 1237 

Longmire C. L. & Rosenbluth M. N., 1956, Phys. Rev., 103, 
507 

Martinez P. V. , Lites, B. W. & Skumanich A., 1997, ApJ, 
474, 810 

Mestel L. & Spitzer L., 1956, MNRAS, 116, 504 (MS56) 
Moll R., Cameron, R. H. & Schiissler M., 2011, A&A, 533, 
A126 

Muthsam, H. J., Kupka, F, Low-Basclli, B., Obertscheider, 

C, Langer, M. & Lenz, P. 2010, New Astron., 15, 46 
Nordlund, A., Stein, R. F., & Asplund, M. 2009, Living 

Rev. Sol. Phys., 6, 2 
Ofman, L. & Thompson, B. J., 2011, ApJ, 734, Lll 
Pandey B. P. & Wardle, M., 2006, astroph/0608008 
Pandey B. P. & Wardle, M., 2008, MNRAS, 385, 2269 
Pandey B. P., Vranjes, J. & Krishan, V., 2008, MNRAS, 
386, 1635 

Pandey B. P. & Wardle, M., 2012a, MNRAS (in press), as- 
troph/arXiv: 1 108.3169vl , 

Pandey B. P. & Wardle, M., 2012b, MNRAS, 423, 222 

Parker, E. N., 1979, Cosmic Magnetic Fields, (Oxford: 
Clarendon Press) 

Parhi,S., Pandey, B. P., Goossens,M., Lakhina G. S., de 
Bruyne P., 1997, Ap&SS, 250, 147 

Parhi,S., Pandey, B. P., Lakhina G. S, Goossens,M. & de 
Bruyne P. 1997, AdSpR, 19, 1891 

Parhi,S., Pandey, B. P., Goossens,M. & Lakhina, G. S. 
1998, IAU, 185, 457 

Petrovic,D., Vranjes,J., Poedts,S., 2007, A&A, 461, 277 

Priest, E. R. 1987, Solar Magnetohydrodynamics, (Dor- 
drecht: D. Reidel) 

Schiissler , M. 1984, A&A, 140, 453 

Schrijver, C. J., Cote/, J., Zwaan C. & Saar, S. H., 1989, 

ApJ, 337, 964 
Schiissler M., 1984, A&A, 140, 453 
Spiegel, E. A., & Veronis, G. 1960, ApJ, 131, 442 
Shelyag, S., Keys, P., Mathioudakis M. & Keenan, F. P., 

2011, A&A, 526, A5 
Solanki, S. K., Keller, C, & Stenflo, J. O. 1987, A&A, 188, 

183 

Soler, R., Terradas, J., Oliver, R., Ballester, J. L. & 

Goossens, M. 2010, ApJ, 712, 875 
Stein, R. F., & Nordlund, A. 1998, ApJ, 499, 914 
Stenflo, J. O., Solanki, S. K., & Harvey, J. W. 1987, A&A, 

171, 305 

Steiner, O. & Rezaei, R., 2012, arXiv:1202.4040vl 
Sykora, J. M., De Pontieu, B. & Hansteen, V., 2012, ApJ, 

753, DOL10.1088/0004-637X/753/2/161 
Vernazza J. E., Avrett E. H. & Loser, R., 1981, ApJS, 45, 

635 

Vranjes, J., Poedts, S. & Pandey, B. P., 2007, Phys. Rev. 

Lett., 98, 049501 
Vranjes, J., Poedts, S., Pandey, B. P. & De Pontieu, B. 

2008, A&A, 478, 553 
Wardle M., 1999, MNRAS, 307, 849 
Wardle M., 2007, Astroph.Sp.Sci., 311, 35 
Wardle M. & Salmeron R., 2011, (as- 

troph/arxiv:1103.3562vl) 



Wedemeyer-Bohm S. & Voort, L. V., 2009, A&A, 507, L9 
Zaqarashvili, T. V., Diaz, A. J., Oliver, R. & Ballester, J. 
L. 2010, A&A, 516, A84. 
Zirker J. B., 1993, Sol. Phys., 147, 47 



APPENDIX A: 

Often, ion-neutral collision frequency is calculated by as- 
suming collision cross-section a = 10~ cm (Khodachenko 
et al 2004). This cross-section is even smaller than H + + H 
collision becomes clear, if we recall that, proton-neutral hy- 
drogen cross-section is well described by the following power 
law at low energies (Glassgold et al. 2005) 



a(E) = g(E 1 ) 



where Ei 



E_y 



(Al) 



0.01 eV and a(E x ) = 1.65 x 
In solar atmosphere, E > Ei 
and p = -1/8. For Model C, VAL81, a « 10~ 14 cm 2 . Thus, 
if we assume that ion-neutral collision is solely due to 
the ionised and neutral hydrogen, the ion-neutral collision 
frequency will be an order of magnitude higher than when 
a « 10~ 15 cm 2 . 

However, lower solar atmosphere is populated mainly 
by the metallic ions and neutral Hydrogen and thus, the 
collision cross-section given by Eq. (Al) underestimates the 
collision frequency. Further, the quantum effects could also 
become important. Thus, it is important to find an ana- 
lytical expression for the momentum transfer cross-section. 
At low energies (< 10 4 K), the elastic cross-section can be 
well approximated by Massey-Mohr cross section (Krstic & 
Schultz 2009) 



omm = 5 x 10 



2/5 



(A2) 



Here Ce is the dipole-dipole coefficient in the interaction po- 
tential which can be expressed in terms of polarisabilities a 
and ionisation energies E of two interacting atoms (Fontana 
1961) 



Ce = 1.5 Qi «2 



Ei E-2 



(A3) 



E\ + E2 

The value of Ce varies between 1.26 x 10~ 60 erg — cm 3 for He 
to 2200 x 10~ 60 erg - cm 3 for Cs (Fontana 1961). Therefore, 
we assume Ce = 100 x 10 _60 erg — cm 3 . Since the cross- 
section is proportional to Cg' 4 , the assumed value of Ce will 
not significantly affect the collision cross-section. The colli- 
sion rate is 

< ctmmv > in = 3.15 x 10 12 v 3/5 cm 3 s _1 (A4) 
which can also be written as 

< ommv > in = 7.91 x 10~ 10 T ' 3 cm 3 s _1 . (A5) 

Here we have assumed v = \JkB T/rrn. The resulting value 
of v in is an order of magnitude smaller than has been as- 
sumed by Vranjes et al. (2008) owing to their arbitrary mass- 
scaling of the cross-section. 

The electron-neutral collision rate is assumed as 



<av > e „= 10 



is ( 128 k B T 
9irm e 



8.28x 10" 10 yTcm 3 s" 



1 -(A6) 
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For similar reasons, the electron-neutral collision frequency 
in the present case is an order of magnitude smaller than 
given by Vranjes et al. (2008); Pandey et al. (2008). 
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